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Abstract 

For two spatially anisotropic, SU(2)-invariant models of frustrated magnets in arbitrary space 
dimension we present a non-perturbative proof of the existence of neutral spin- 1/2 excitations 
(spinons). In one model the frustration is static and based on fine tuning of the coupling constants, 
whereas in the other it is dynamic and does not require adjusting of the model parameters. For both 
models we derive a low-energy effective action which does not contain any constraints. Though our 
models admit the standard gauge theory treatment, we follow an alternative approach based on 
Abelian and non-Abelian bosonization. We prove the existence of propagating spin-1/2 excitations 
(spinons) and consider in detail certain exactly solvable limits. A qualitative discussion of the most 
general case is also presented. 

PACS numbers: PACS No: 71.10.Pm, 72.80.Sk 
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I. INTRODUCTION 



Excitations with fractional quantum numbers in general and neutral excitations with 
spin-1/2 in particular constitute an essential ingredient of many theoretical approaches 
aimed to explain the non-Fermi liquid behavior of quasi-two-dimensional copper oxide ma- 
terials. These approaches try to generalize for higher dimensions such mechanisms as spin- 
charge separation and quantum number fractionalization, the phenomena well understood 
and adequately described for one- dimensional systems (see recent review articles (T| , , || 
and references therein). 

As a part of the general program to construct a theory of copper oxides, there have been 
many attempts to find neutral spin-1/2 excitations (spinons) in purely magnetic systems, 
that is to find higher-dimensional analogues of ID spin S=l/2 Heisenberg antiferromagnet. 
Such generalization was qualitatively outlined by Anderson J| in the form of the famous 
Resonant Valence Bond (RVB) state. This is a spin-liquid state which breaks neither trans- 
lational nor spin rotational symmetry. Building on this proposal, Kivelson, Rokhsar and 
Sethna showed |5] that such state, if exists, must support neutral spin-1/2 excitations. 

Despite more than a decade of strong efforts, no realization of an RVB state supporting 
excitations with fractional spin in D > 1 Heisenberg magnets with short-range exchange 
interactions has been found in models with interactions of the Heisenberg type. We exclude 
from the consideration models of quantum dimers, where such excitations have been shown 
to exist ||, and other recently suggested models with fractionalized excitations, which are 
not microscopic electronic models (see, for example as well as generalizations for 

symmetries higher than SU(2). We also do not consider magnets with incommensurate 
ordering (such as a spiral one) which also support spin-1/2 excitations P|,[^0[|. Strong 
qualitative arguments against the RVB scenario have been presented by Read and Sachdev 
PH , wno argued that the most likely mechanism for disordering an antiferromagnetic state 



is the spin-Peierls one. Here, although spin rotational symmetry is preserved, translational 
symmetry is broken since the spin system undergoes an explicit dimerization at a finite 
temperature. That makes the existence of deconfined spinons in the low-temperature phase 
impossible. 

An interesting by-product of the RVB scenario is Strong Coupling Compact Quantum 
Electrodynamics (cQED) field theory. This gauge theory describes the so-called 7r-flux RVB 
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state in the continuum limit |12| (see also [|TI|,JT4|]), under the assumption that the system 



does not develop any spin-Peierls ordering. The gauge theory has infinite bare coupling and 
therefore is difficult to study. All existing analysis relies on the results obtained for the SU(N) 
or Sp(2N) generalizations of this theory based on the 1/N-expansion, as in |TT[],|12||,[|15|],||16|. 



In this paper we describe two D-dimensional (D > 1) models which exhibit neutral spin- 
1/2 excitations (spinons) and apparently represent realizations of the 7r-flux RVB state. 
The paper is organized as follows. In Section 2 we introduce the models. They are spatially 
anisotropic and actually represent collections of weakly coupled chains. In all of these 
models, the naively strongest interactions between the chains, namely those that couple 
the antiferromagnetic (with momentum close to 7r) fluctuations on neighbouring chains, are 
frustrated and so can be set to zero in the Hamiltonian. In this case, the physics is governed 
by a marginally relevant interaction between the long-wavelength magnetic fluctuations. In 
Section 3 we discuss results we managed to obtain for infinite number of chains. These can 
be summarized as follows: 

• we give a proof of the existence of spin-1/2 excitations in the limit of infinite number 
of chains; 

• we demonstrate that in the low-energy limit the effective spin Hamiltonian decouples 
into two commuting parts describing sectors with different parity; 

• we show that the models are stable against an antiferromagnetic phase transition; 

• for the models in question we demonstrate existence of a zero temperature phase 
transition (spontaneous transverse dimerization) . 

• we show that the ground state is degenerate; for periodic transverse boundary condi- 
tions the degeneracy is equal to 2 2N ~ 1 where 2N is the number of chains. 

To achieve a more detailed understanding of the model, its spectrum and correlation 
functions, we study some limiting cases. Thus in Section 4 we discuss two exactly solvable 
cases: two and four chains with periodic boundary conditions in the transverse direction. In 
Section 5 we consider another solvable case: the model of three chains with open boundary 
conditions. In all these cases magnetic excitations carry spin 1/2, in agreement with the 
general result. Systems with more than two chains also display non-magnetic (singlet) modes 
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and a topological order. We show that this order is similar to the one which exists in the 
ID Quantum Ising model. 

Further, in Section 5 we introduce a two-dimensional model with a hierarchy of interac- 
tions. We use this, apparently artificial but solvable model as a substitute for the original 
(spatially uniform) one with a hope to extract more detailed information about the structure 
of the Hilbert space. The hierarchical model is solved by a cluster expansion. The results 
indicate that about half of the Hilbert space is occupied by singlet modes. The ground state 
remains massively degenerate with the ground state entropy proportional to the number of 
chains. 

Having in view future numerical simulations for our models, we discuss in Section 6 
finite-size effects for systems of a finite number of chains. In Section 7 we use the acquired 
knowledge to conjecture the properties of correlation functions in the thermodynamic limit. 
Our conclusions are summarized in Section 8. 

II. THE MODELS 

As in the original suggestion by Anderson [H, the element most essential for our construc- 
tion is frustration. Another ingredient of our approach is strong spatial anisotropy: so far 
we have only been able to deal with models where the exchange interaction in one direction 
is much greater than in the others. Thus it is proper to describe our models as assemblies of 
weakly coupled chains. Our results also suggest an alternative approach to Strong Coupling 
cQED. We will not elaborate on this analogy, though, postponing this interesting question 
for future studies. 

The first model is a spin-1/2 Heisenberg magnet on an anisotropic Confederate Flag (CF) 
lattice. In what follows, we will be working with a two-dimensional version of this model. 
The interaction pattern depicted on Fig. la; one can easily generalize this construction to 
three dimensions. The CF model Hamiltonian is given by 

Hacb = H WI|S 

j,n { fi=±l ) 

where S J>n are spin-1/2 operators, and Jy, J r , > 0. Since the interchain couplings (J r , Jj) 
are much smaller than the exchange along the chains (Jii) it is legitimate to adopt a con- 
tinuum description of individual chains. In this description, the local spin densities are 
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represented as sums of the smooth and staggered parts: 



Sj,n/ao -> S n (x) = M n (x) + (-l) J N n (x), x = ja Q , (2) 

do being the lattice spacing in the chain direction. The low-energy dynamics of the spin- 1/2 
Heisenberg antiferromagnet 

H id = J\\ XX S i S i+i) ( 3 ) 

3 

is described by the SUi(2) Wess-Zumino-Novikov-Witten model. The latter Hamiltonian 
can be written in terms of the so-called chiral vector current operators, J and J, satisfying 
the level k = 1 Kac-Moody algebra (this approach has been described in a vast number of 
publications; see for a review [|T7|] or fT5|): 

~:(J.J): + :(J.J):1+..., (4) 



Hid -> / d.r 



where v = irJ\\ao/2 is the spin velocity and the dots stand for a marginally irrelevant 
perturbation which will be discarded in what follows. It is remarkable that the smooth part 
of magnetization, 

M = J + J, (5) 

and the spin current, 

j = «(J-J), (6) 

are locally expressed in terms of the chiral currents. 

In the CF model, the exchange is frustrated in the direction perpendicular to the chains 
and can be fine tuned to make its Fourier transform vanish at q» = n. This is achieved when 
the rung and plaquette-diagonal coupling constants satisfy the relation J r = 2Ja, in which 
case the direct interchain interaction between staggered magnetizations, N n (x) ■ N n+1 (x), 
is completely eliminated. The absence of this, strongly relevant perturbation is the most 
important property of our model. The triangular lattice is excluded from the consideration. 
It was demonstrated in |19) that the coupling of the staggered parts of magnetization cannot 



be removed completely in that case: parity-breaking ( "twist" ) terms surviving the continuum 
limit make the analysis very complicated. 

Apart from the collective spin excitations, the second model also involves charge degrees 
of freedom. It is a quasi-one- dimensional Kondo- Heisenberg model (see Fig. lb) which has 



a) 2J g = J r « J b) 

FIG. 1: a.) Exchange interactions pattern for CF model, b.) Lattice of stripes for the Kondo- 
Heisenberg model. 



been already discussed in the context of theory of stripes P0fl,pi| ( see a l so 0] where the 
two chain case is discussed). It is assumed that neighboring chains have different band 
filling such that metallic chains are surrounded by insulating spin-1/2 magnetic chains. 
In order to eliminate the electron tunneling across the magnetic chain, one needs to have 
magnetic stripes consisting of at least three chains. At low energies one can consider such 
stripe as an effective single spin-1/2 chain. Here the frustration is dynamical: the electrons 
living on metallic chains do not experience back-scattering on magnetic excitations due to 
incommensurability of the Fermi wave vectors: 2kp ^ ir. It is easy to demonstrate (see, for 



instance |20||) that the charge excitations in this model decouple and remain one-dimensional 
(this is true at least in the first approximation when one discards various virtual processes). 
On the other hand, the spin sector is described by the same effective Hamiltonian as in the 
CF model. 

Using the continuum description of individual chains in the interchain exchange inter- 
action, based on the asymptotic representation (Q) of the spin operators, we arrive at the 
following effective Hamiltonian: 

2N 



n=l 



HlD,n + — X! + J) n ' + 
1 n=±l 



(7) 



where Hid is given by Eq.(f|). Here 7 is determined by the lattice Hamiltonian; in the CF 
model with interactions chosen as in Fig. la we have 7 = 2J r a . For the Kondo-Heisenberg 
model shown in Fig. lb, Eq.(^) is modified in such a way that even (metallic) and odd chains 
(effective spin-1/2 chains representing a stripe of three coupled Heisenberg chains) have 
different velocities. This detail, however, does not lead to any qualitative changes. 

The model (0) is closely related to the gauge theories extensively studied in the context 



of strongly correlated systems. This analogy is briefly discussed in Appendix A. In the rest 
of this paper, however, we will follow a different route. Namely, we will return to Eqs.(^), 
(0) and employ the Abelian bosonization which proves particularly useful in revealing the 
topological nature of elementary excitations and determining their quantum numbers. 

III. PROOF OF EXISTENCE OF THE SPIN-1/2 EXCITATIONS. INFINITE 
NUMBER OF CHAINS 



It is well known (see e.g. |Tj|) that the chiral SUi(2) currents can be faithfully represented 
in terms of holomorphic (left) and antiholomorphic (right) scalar fields, <p and tp: 

J n = 4=^n, Jn = (8) 

Jl = -4=^». Jn = -t-e*^. (9) 
v27r 2na 



Here d, d = \{v l d T ^\d x ), a is the ultraviolet cut-off of the bosonic theory, and the fields (p 
and <p are governed by the chiral Gaussian actions f2"3fl, ||24| (the latter work used this form 



of bosonization in the context of the theory of edge states in Quantum Hall Effect). The 
bosonized low-energy effective action can then be suitably represented as 

N 



S=j: [ drdx [£+ + £;+ £™\, (10) 

n=l J 

£+ = d x (p2n(d x ip2n ~ W~ l d T <p 2 n) + d x Lp 2n+ i(d x Lp 2n+ i + IV' 1 d T (p 2n +l) + 

+ 7 E ((27ra)~ 2 cos[v / 8vr(<^2n + <^2n+ M )] + ^-d x <p 2n d x <p 2n+ A , (11) 
M =±i L 2vr J 

C~ = - 0), (12) 

£> % n =7E (2vr«) _2 cos[v / 87r(v5 2 n. - <^2n+ M )] + —d x <p 2n d x ip 2n+fl \ + (<p -> <p). (13) 
M =±i 1 2vr J 

This is the form of the action in which the requirement of vanishing charge current (J** = 0), 
imposed in the gauge field RVB approach, is explicitly resolved. 

According to the definitions (||), @, the total spin projection S z is equal to 

1 r°° 

S Z = ^=J2 dx(d xVn + d x <p n ). (14) 

V 27T n J -co 

This is a general definition which does not assume that the fields <p and (p are independent. 
In fact , the interchain interaction that flows to strong coupling in the low-energy limit freezes 
certain combinations of the fields with opposite chiralities and thus makes them coupled. 



In the first loops the Renormalization Group (RG) equations for the model of infinite 
number of chains coincide with the equations for two chains. At 7 > the interaction of 
currents with different chirality in (|11|), fll^D is marginally relevant and reaches the strong 
coupling at an energy scale 

A = C (v 7 ) 1/2 exp( -7^/7), (15) 

where C is a number. At the same time, the (Lorentz-noninvariant) perturbation C tnt , which 
is responsible for the velocity renormalization, is irrelevant and flows only in the presence 
of the other interaction. Thus it is likely that this interaction remains weak even when the 
relevant one reaches the strong coupling. For this reason, in the situation when the bare 
constant 7 <C 1, we deem it possible to neglect C mt . This results in an approximate Lorentz 
invariance in (l+l)-dimensional space-time which plays an important role in the analysis 
that follows. Without C mt , the action splits into two independent, chirally asymmetric 
sectors, S = S + + S~. The two sectors are mapped onto each other under the reversal of 
the chiralities, ^ n <-> (^„, or a shift by one lattice spacing in the transverse direction. Hence 
the total symmetry becomes [SU(2)]+(g>[SU(2)]_ and the excitations in this model carry two 
quantum numbers: spin and parity. 

Notice that the multi-chain effective action (|T0|) features a local Z2 symmetry already 



mentioned for the two-chain problem p5| : as opposed to the original lattice model which 
displays global translational invariance, the effective Hamiltonian of the continuum theory 
(0) is fully expressed in terms of the vector currents and, therefore, is invariant under 
independent translations along the chains by one lattice spacing. These translations are 
generated by the shifts 



<£> n -> <£ n ± W -, ip n ^lf n , (16) 

or 

V?n -> <fn, <fn ~» ± J^- (17) 

The development of strong-coupling regime in the low-energy limit leads to a spontaneous 
breakdown of this symmetry. Indeed, in the ground state the field combinations <p n + {p n +n 
get frozen at the values (see Appendix B for the details) 

^» + ft*.M = ^f+yf™S°. ™^=0,±1,±2,... (18) 
For an even number of chains and with periodic boundary conditions in the transverse 
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direction the integers m^' satisfy 

= 5>(->. (19) 

n n 

The sets of integers {m^'} subject to the condition ( |T9"D label different degenerate ground 
states, each of them breaking the local translational symmetry (|TB|), (|T7|). The existence 
of stable degenerate minima of the potential implies the existence of massive topological 
excitations, solitons, interpolating between adjacent vacua. Using (|l4l) and ( |l8|) one can 
represent the spin of a given state as 

s z = ^K +) N- m i +) M]. (20) 

Thus, the minimal nonzero value of the spin is 1/2. This is the spin carried by the topological 
solitons. 

We emphasise that this statement is based on the exact symmetries of multi-chain Hamil- 
tonian fllCf ) and does not depend on the adopted approximations. Within the approximation 
that neglects C mt there are solitons and antisolitons with a well defined (+) and (-) parity. 

We would like to briefly comment on the above results. The factor v^vr that appears in the 
arguments of the cosines in Eqs.([TT|), (|12|), representing a marginal current-current interchain 
interaction, is crucial. Recalling the bosonization rules for the staggered magnetization of a 
single spin-1/2 chain (see e.g. |T8|) 



N n ~ (cos v / 2tt(^„ - <p n ), smV2n(<p n - tp n ), - sin \/2ir(ip n + tp n ), ) (21) 

we realize that interchain coupling between staggered magnetizations would give rise a 
strongly relevant perturbation containing cosines with \[2tx. If such terms were present 
in the action, the period of the potential and, hence, the topological quantum numbers of 
the excitations, would be doubled as compared to our case, and thus the spinons would be 
confined. However, by construction these terms do not appear in our models, and this is 
why the interchain interactions do not lead to spinon confinement. This argument does not 
require neglecting L int because the latter term also contains only y/87t cosines. 

Thus we are confident that in the models in question spin-1/2 excitations do exist. We 
also may rest assured that these excitations do not remain confined to individual chains, as 
it happens, for instance in the so-called 'sliding' Luttinger liquid phase in arrays of crossed 
chains ||27||. In that latter case, the spinons survive due to irrelevance of the interchain 



interactions. Contrary to that scenario, in our case the interchain coulings remain relevant 
making the system truly multi-dimensional. 



A. A failure of the 'most obvious' solution 

The following calculation suggests that freezing of the fields does not imply that one 
can expand around the vacuum configurations. Naively, one would adopt a self-consistent 
harmonic approximation (SCHA) in which the cosine term is expanded around the minimum: 

1 



cos[v / 8vr(^ n + v?n+ M )] « const + ^p[<%„ + <Pn +A1 )] 2 '. (22) 

Here 5 denotes a deviation from the vacuum value. SCHA assumes that the transverse 
stiffness p should be determined self-consistently: 

p = — y{exp[W8n((p n + (pn+n)]}- (23) 

We shall demonstrate, however, that such procedure yields zero p. 

The effective action for fluctuations around the minima is (here we set v — 1) 

/ 



(24) 



q(q + ko) + p, p cos k \ I a p 
I pcosk, q(q-iu)+pj I (3 P 

where a = Sip, (3 = 5ip, p = (u, q,k), q and k being the longitudinal and transverse momenta, 
respectively. Self-consistency condition (|23|) becomes 

q 2 + p(l — cos k) 
0J 2 q 2 + (q 2 + p) 2 — p 2 cos 2 k 

(25) 



lnp ~ -[G aa + Gpp + 2G Q/3 ](r = 0,x = 0) ~ J dtodqdk 
dkdq q 2 + p(l — cos k) 



(q 2 + p) 2 — p 2 cos 2 k 

Since the latter integral diverges, the suggested naive solution is not self-consistent. Its fail- 
ure originates from the chiral structure of the Gaussian part of action ( |i~T|) which includes 
the first power of the time gradient. For the conventional (non-chiral) Gaussian model a 
similar solution would be self-consistent. One immediate consequence of it would be the ap- 
pearance of an average staggered magnetization. Thus the present calculation demonstrates 
stability against antifferomagnetism. 
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B. Order parameters 



As we have already seen, the fact that the combinations of the fields (f2 n + y?2n+^ and 
V?2n + V?2n+^ get locked at + n, {n <G Z) implies that the ground state of the system 
is characterized by a spontaneously broken local translational invariance (|16j), (|17D. The 
most obvious order parameters (OP) in the two sectors of the model (|T0|), whose structure 
directly follows from the form of the interaction, are given by 

7^ {+) = i(exp[iA/27r(v9 2n + (p 2n +^)]), 

= i(exp[i V / 2^(^2n + <^2n+ M )]> (26) 

The above OPs are real (see Eq. ([T8|) ), defined on the links of the lattice and are non-local 
in terms of original spins. Moreover, since the OPs include chiral fields, they are expressed 
(nonlocally!) in terms of magnetization density and spin currents. As is shown in Appendix 
B, for different pairs of neighboring chains, the signs of the order parameters are uncorrelated. 
This leads to a massive gound state degeneracy 2 2N ~ 1 (recall that 2N is the total number 
of chains). 

The non-locality of the order parameter suggests a possibility of topological order. This 



supports the original idea by Wiegmann that topological order exists in RVB states [p9fl ,fl30 
Characteristic features of such order is non-locality of the order parameter, ground state de- 
generacy (in a closed system) and the existence of edge states (in an open system). The 
ground state degeneracy depends on the topology of the manifold (hence the name 'topo- 
logical'). A relative stability of the state with respect to perturbations is protected by the 
spectral gap. One well known example of topologically ordered states is provided by incom- 



pressible Fractional Quantum Hall states (see the discussion in [31] and references therein). 
In our case the degeneracy of the ground state is exponential in the number of chains. This 
fact alone indicates that these models have a topological order of a different type. This 
order is closely related to the one existing in 2D Ising model, multi-channel Kondo models 
and models with scattering matrices of the RSOS-type. 

The first examples of non-local order parameters are given by the disorder field of 2D Ising 
model |52| and the dual field exponentials in 2D XY model|l3|,|28|]. These fields acquire 



non-zero expectation values in disordered phases of the corresponding models. According 



to P^|,P5|,P^|, [P7P,|]I5|, such order parameters (also called 'string' order parameters) exist 
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in various massive phases of ID spin liquids. All these liquids are related to either spin 
S=l Heisenberg chain or two-leg spin 1/2 Heisenberg ladders. As was demonstrated in 
58||,[P7||)[P9| (see also jn§), after the Wick rotation these systems can be described as a 



collection of 2D Ising models. The evidence that the order is really topological was presented 
by Affleck [[KJ, who proved that ID spin liquids have edge, or rather end, states: open chains 
have free spin-1/2 states located at the chain's ends. 

We postpone a more detailed discussion of the issue of topological nature of the order till 
the further Sections. 

Order parameters (|26D change their sign on soliton configurations and, therefore, vanish 
at finite temperatures when the densities of solitons and antisolitons are finite and equal. 
This mechanism is a sufficient (though not a necessary) one to destroy the order parameter 
at finite T. The sign changes may also originate from fluctuations in non-magnetic channels. 
Later we will discuss an exactly solvable example of four coupled chains in which the latter 
scenario is realized. The temperature dependence of the correlation length depends on 
what mechanism is realized. If the order parameter is destroyed only by (anti) solitons, the 
correlation length is exponentially large in 1/T. This is because the topological excitations 
have a finite spectral gap and therefore their density is exponentially small: n ~ exp(— A/T). 
In this case, the correlation length is defined as the average spacing between the solitons: 
£ ~ rT x l D (D is the space dimensionality). 

As we shall see, the singlet excitations either have a very small gap or even no gap at 
all (we do not know whether the latter is true). Then the order parameter is governed 
by another correlation length quite distinct from the above. The density of non-magnetic 
excitations above their gap depends on T as a power law. Such situation would correspond 
to a quantum critical point with a power law correlation length. 

Turning to physical fields, we note that the products of the staggered magnetizations on 
neighboring chains representing the transverse dimerization order parameters (see Appendix 
B) are expressed as bilinear combinations of the true order parameters: 

(NLNL+,) = ^H C0S {V2^2n + ^2n +M + ^2n + </W}> = j[T«T« + c.c],(27) 

= (N y 2n N y 2n+ ») 
■(exp[iv / 27r(^2n + <^2n+ M )] exp[-iv / 27r(<^2n + <P2n+»))} 



8na 

- A {^ +] [T^r + c.c.}. (28) 
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Since T* = T we have 

implying that the SU(2) symmetry remains unbroken in the ground state. 

Now let us consider the magnetic correlation functions. The order parameter field of 
ID spin-1/2 chain is not just a staggered magnetization; it is a 2x2 SU(2) matrix g whose 
entries include not only vector components of the magnetization n a (a = x,y,z), but also 
the dimerization field e = (— l) n (S Tt S n+ i): 

g(x) = e(x)I + ia a N a (x). (29) 

It can also be represented in a factorized bosonic form: 

<W = ~^ C ™' : expf-iv^o-y? + a'Cp)} := C atJ >z c z a >, 

C aa , = e Mi-™')/\ (30) 

where 

z a = exp[iay/2Trtp], z a = exp[— iaV2n(p], (a = ±1). (31) 

As follows from the factorization of the low-energy effective action, if we consider correlation 
functions on chains with the same parity (for instance, both even) z correlate only with z's 
and z with z. For Green's functions on chains with different parity z correlate with z and 
not with z. Therefore the correlation functions of the staggered magnetizations and energy 
densities factorize: 

(g ab {r, x- 2n)g+ (0, 0; 2m)) = (z a (r, x; 2n)z*(0, 0; 2m)){z b (r, x; 2n)z* c (0, 0; 2m)) 
= 5 ad 5 bc G nm (T, x)G nm (r, -x), (32) 
(g ab (r, x- 2n)g+(0, 0; 2m + 1)) = (z a (r, x; 2n)z* c {0, 0; 2m))(z b (r, x; 2n)z* d (0, 0; 2m + 1)) 
= 8 ac 5 bd D nm (T,x)D nm (-T,x). (33) 

Various components of the z-s are conveniently assembled into vectors and co-vectors: 



o-,2n' Z cr,2n+l I ) 



( Z \ 



(34) 
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for the (+) sector and 

( Z l,2n+n K,2n) J ! I (35) 




for the (-) sector. Then the corresponding correlation functions become matrix elements of 
the 2x2 Green's function: 

/ 



G 



(36) 



G(uj,k z ;k ± ) D(u,k z ;k ± ) 
\ D(-u, k z \ k±) G(aj,-k z ,k x ) 
An important fact is that both the function D and fcj_-dependence of G are non-perturbative 
effects: they do not appear in any order of perturbation theory in 7 and, as we shall see 
later, are proportional to exp(— ttv /j). 

IV. EXACTLY SOLVABLE CASES: TWO AND FOUR CHAINS. 

From now on we will discuss the (+) and (— ) sectors independently. To distinguish 
between the amplitude of the RVB order parameter and various mass gaps, we shall denote 
spectral gaps by the letter m. 

To check the ideas outlined in the previous Section, we shall discuss some exactly solvable 
cases. Consider a two-dimensional array of an even number (2N) of chains with periodic 
boundary conditions. In this case we can arrange chains in pairs. Exact solution is possible 
in two cases: N = 1 and N =2. 

Throughout this Section we shall use the so-called form-factor approach to calculate 
leading asymptotics of the correlation functions. This method has been widely applied to 
massive integrable field theories and there are excellent review articles |^2[ and even books 



about it. The form-factor approach uses information extracted from the scattering the- 
ory and symmetry properties of physical operators to calculate matrix elements of these 
operators between various exact eigenstates. To estimate the correlation functions one uses 
the Lehmann expansion. This approach is well defined for theories with massive excita- 
tions. In one dimension, for reasons of kinematics, matrix elements of multi-particles states 
become smaller and smaller when the number of excited particles increases. This feature 
greatly improves the convergence of the Lehmann expansion and allows one to obtain an ac- 
curate description of the correlation functions with only a few (sometimes even one!) matrix 
elements involved. 
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A. N =1 (two chains) 



Though the number of chains is really small, this case contains certain interesting features 
which may shed light on more general cases. 

For N =1 each sector with a given parity is equivalent to the isospin sector of the SU(2) 



invariant Thirring model, as it was shown in [25]. The excitations are massive solitons with 



spin 1/2. There are no other modes (singlet ones), as it must be for the case of large N. 
For this reason this case is not representative, as we have already said. The massive modes 
correspond to amplitude fluctuations of A. 

The correlation functions of chiral bosonic exponents for the Thirring model can be ex- 



tracted from the recent paper by Lukyanov and Zamolodchikov |g4j (see also |4l|j). In par 



ticular, for the correlation functions of spinon fields ( ftq) and for the correlators of staggered 
magnetizations we have 

G(r,x) = (ei^^e-iVSFvCo.o)^ „ 1 e xp(-mpi 2 ) 

y/iVT — X 

^VSF* 1 (r, iB ) e -i^ 1 (o,o)v „ ± exp( _ 2mpi2 ) 
D(t,x) = {j^™j^W)~ m WK (mp 12 ) 

^ e iv^i(T,x) e iv^ 2 (o,o)^ = ZQmK 2( mpri ) ( 37) 



where p = y t 2 + (x/v) 2 and the mass gap m = A given by Eq. ([IS]) . In writing down these 
formulae we used the single mode approximation for the spinon operators (for the details 
we again refer the reader to [pBf): 

z+(t, x) = Zl /2 m 1/4 J d6e e/i e- ixms ' mhe [e- mTCOshe Z+{8) + e mTCOshd Z_ a {9)\ (38) 
z a (r, x) = Zl /2 m l/i J d8e- e/i e- ixmsinhe [e- mTCOshd Z+(6) + e mTCOshe Z_ a {6)] (39) 

where Z±(Z±) are annihilation (creation) operators of solitons and antisolitons and Zq is 
a numerical constant. Approximation ( |3"9"D fails at small distances; inclusion of terms with 
multiple soliton production removes the logarithmic singularity of K . 



The two-point correlation function of currents were calculated in [45 . 

Eq . (|37D supports the earlier claim made at the end of the previous section, that the 
interchain correlation function of the staggered magnetization is non-analytic in 7 and thus 
vanishes in all orders of perturbation theory. 
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B. N =2 (four chains) 



In the N =2 case the interaction can be written as 

(Ji + J 3 )(J 2 + J 4 ) + (J 2 + J 4 )(Ji + Js) (40) 
A sum of two k = 1 SU(2) currents is the k = 2 current; moreover, according to HB 



the sum of two SUi(2) WZNW models (the central charge 2) can be represented as the 
SU2(2) WZNW model with central charge 3/2 and plus one massless Majorana fermion (a 
critical Ising model) with central charge 1/2. Using the results of |||| we rewrite the entire 
Hamiltonian density (fD as follows (here only the (+)-parity part is written): 

1~~("massless Umassive 

iv iv 

T-tmassless = —^XO^xXO + jXod x XO (41) 
7T1) — — — M) 

H massive = — (: I • I : + : I • I :) + 7 I • I = -(- X a d x X a + X°^X°) - \(x a X a ? (42) 
where a = 1, 2, 3 and 

I = Ji + J 3 , I = J 2 + J 4 

Ji,3 = ^{±XoX+ \{X x X]}, J2,4 = ^ |±XoX+ 2^ x X]} (43) 

The fields x, X stand for real (Majorana) fermions. 

Eq.flilD describes a critical Ising model; the corresponding excitations are gapless and 
non-magnetic; they appear in the sectors with both parities. The criticality is an artefact 
of the periodic boundary conditions in the transverse direction. 

The massive sector can be described using three equivalent representations. One of those 
is the SU2(2) WZNW model perturbed by a current-current interaction (the ultraviolet 
central charge C\jy = 3/2). Here the fundamental field is the SU(2) matrix g. Another 
representation is the 0(3) Gross- Neveu (GN) model where fundamental fields are Majorana 
fermions (0). The third representation is three critical Ising models coupled by the energy 
density operators ef. 

3 

3 = 3lsing{Vi) ~^Y.J drdx € i e J ( 44 ) 
i=l i>j 

In the latter case the fundamental fields are order parameter fields of the Ising models or, 
via the duality relationships, the corresponding disorder parameter fields. 
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An external magnetic field couples to the sum of all currents and therefore only to the 
fermions of the 0(3) GN model: 

H magn = ih • ([£ X ^ + [| X *]) (45) 

Since this model has a spectral gap, the net magnetization does not appear until the mag- 
nitude of the field reaches the value of the gap. 

The Thermodynamic Bethe Anzats equations for 0(3) GN model were obtained by one 
of the authors HTJ; some of the correlation functions were calculated by Smirnov A 
very illuminating discussion of the symmetries and the excitation spectrum of this model 
can be found in the recent paper by Fendley and Saleur ED . 



The massive magnetic excitations are particles with a non-Abelian statistics, as in the 
Pfaffian state of Quantum Hall effect fl5"0| . An unusual fact about the 0(3) GN model is 



that its ground state is doubly degenerate. This degeneracy is of topological nature being 
related to properties of zero modes formed on solitons. The picture simplifies drastically if 
one takes into account that 0(3) GN model is equivalent to Supersymmetric sine-Gordon 
model (SUSY SG) 



£ = \(d^) 2 + \xi^x + M 2 sin 2 (/5$) + iMxxcos(/3$) (46) 

with (3 2 = 4tt. For SUSY SG model there is a semiclassical limit (3 2 <C 1, where things are 
easier to grasp. In this limit the spectrum is determined by the minima of the bosonic part 
of the potential. There are two sets of minima: ones where cos(/3$) = 1 and ones where it 
is -1. This sign difference does not affect the energy spectrum of the Majorana fermion but 
affects expectation values of the Ising model order and disorder parameters. Therefore the 
vacua with different sign of cos(/3$) are physically distinguishable. It is clear that the fact 
of this degeneracy is related to (i) odd number of Majorana fermions species present and (ii) 
the existence of solitons. Both these features are very general and are likely to survive in 
a finite parameter domain. The described situation should be opposed to what happens in 
the conventional sine-Gordon model, where the potential is just cos(/3$) and the degeneracy 
between its minima is lifted in the ground state. 

There are two descriptions of the spectrum. The first one was suggested by Zamolodchikov 



51]. This description abandons the notion of asymptotic particle states and, therefore, leads 



to difficulties when being applied to the correlation functions. The second approach was 
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suggested in |52| and used in |48| to calculate the correlation functions. In this approach 



solitons are treated as particles carrying two quantum numbers - an SU(2) spin a =f, j and 
an isotopic number p = ±. The isotopic part of the K-soliton Hilbert space is truncated, that 
is not all combinatorically possible 2 K states are allowed. The number of allowed states is 
2[ ftr / 2 l; in fact, the isotopic space of K solitons is isomorphic to the irreducible representation 
of ^-dimensional Clifford algebra (in other words, one may say that each soliton carries a 7 
matrix - a zero mode of the Majorana fermion). In particular, the two-soliton wave function 
must be an isotopic singlet. 

This approach allows one to use the conventional scattering matrix description at the end 
projecting out unwanted states. The two-particle S-matrix is given by the tensor product 

S(9 12 ) = S IT m(0u) ® S^(6 12 ) (47) 

where the first S-matrix corresponds to the SU(2) invariant (isotropic) Thirring model and 
second one is the S-matrix of solitons of the sine-Gordon model with a particular anisotropy 
£ = 47r (in Smirnov's notations £ = 7r/5 2 /(87r — /3 2 )). 

The order parameters are local operators both in the Ising and the WZNW representa- 
tions: 



T 12 = exp[i v / 27(^i + <p 2 )\ = e^+e 1 ^*- 

= (000-3 + i/ioA i 3)(o"i0'2 + iAtijU 2 ) = 0"o Tn) + i/i Tr<? + 
T14 = exp[iv / 2T(v?i + £4)] = e^+e 1 ^®- = (a /i 3 + {^2 + inm) (48) 

where <7j,/ij are order and disorder parameters of the corresponding Ising models. The 
spontaneous mass generation freezes <Tj, (i = 1,2,3), but leaves <7o,/Uo critical. Therefore 
the most singular parts of the above order parameters 

T12 ~ a , T M ~ yUo (49) 

still have power law correlations at T = 0. The correlation length atT^0is£~l/T and 
is not determined by the solitons - a possibility we have mentioned in Section IIIB. 

The chiral exponents entering in the expressions for the stagerred magnetization are 
non-local: 



expfiv^v^i] = C ab a a {z)g b (z, z) 
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(50) 



where a a are chiral vertex operators of the critical Ising model (see |53 for details). The 
operators g in the ultraviolet limit become vertex operators of the S=l/2 tensor operators 
of k = 2 WZNW model. 

In order to understand the subsequent calculations of the correlation functions the reader 
have to keep in mind two principal facts about excitations in integrable models. The fact 
number one is that these excitations cannot be unambiguously classified as fermions, bosons, 
semions etc.. Their commutation relations are determined by their S-matrix and therefore 



depend on the momenta of the particles (see, for example[56|). The fact number two is that 
the creation and annihilation operators of elementary excitations in integrable theories are 
usually strongly non-local in terms of the bare creation and annihilation operators. Hence 
these excitations are extended objects and therefore do not belong to any particular represen- 
tation of the Lorentz group. Therefore their Lorentz spin is not fixed and the same operator 
can represent physical fields with different Lorentz spin depending on circumstances. 

Correlation functions of the currents 

According to Eq.(^) there are two parts in a given current operator, J_ and I. The 
first one containing a product of the gapless fermion xo an d the 0(3) GN model fermion is 
not conserved. A GN fermion is a convolution of two solitons, each of them entering with 
Lorentz spin 1/4 to make the total spin 1/2. Superficially it may appear that the second 
term (the conserved current) is a convolution of four solitons, but this is incorrect: the 
minimal matrix element contains two solitons which enter here with Lorentz spin 1/2. Since 
Xo remains massless, the leading asymptotics of the correlation function of non-conserved 
currents (J~J~) with the threshold at s — 2m (s 2 = uj 2 — (vq) 2 ) exists only for the currents 
with the same chirality (here we put v — 1): 

(r(r,x)J b _(0,0)} =5 ab (r + ix)- 1 x (51) 
J d0id0 2 |F 1 (0 12 )| 2 e (ei+02)/2 exp[-mr(cosh^i + cosh# 2 ) + imar(sinh0i + sinh0 2 )] 

which gives 

( . r ~ lg) / deme^K^mpcoshO) (52) 

(r + ix) v t z + x z J 

where F\ can be calculated. The Fourier transform 

%m(J-J-) M = f^±l) L(s 2 ), *m{J-J-) M = (^) L(s 2 ), 
\uj - q J \u + qj 
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/(\f) 
r Ms 2 -4m 2 cosh 2 B)\F X (B) I 2 (53) 
cosh 6 1 

where 9h{x) is the Heaviside function. 

Now let us consider the pair correlation function of the conserved 0(3) currents / = 

J\ + J 3 , 7 = 7 2 + 7 4 . Their matrix elements were calculated in j48|. In the scheme adopted 

in that paper a soliton carries two quantum numbers - an SU(2) spin a ="[*, { and an isotopic 

number p — ±. The isotopic part of the multisoliton Hilbert space is truncated such that 

the two soliton wave function must be an isotopic singlet. Due to the SU(2) symmetry, it 

is sufficient to have an expression for one component of the current, for example I 3 . The 

matrix element into a state of two solitons with rapidities 61,62 and quantum numbers 

(er,p)i, (ct,p) 2 is given by 

la ( \ a I Mr*ln\ a cot ^^/^)Coo{0i 2 )Ci{0i 2 ) 
^{a uPl )-6 2 ^ P2 )\I [0> = de« ( g i2 + i7r)cosh[ | ( g i2 + i7r)] >< 

[|U> + Ut>]® [I + ->-!-+>] (54) 

where 

C*(0) = sinh(0/2)G(0), 
Cfc(6>) = exp 

and 



00 sin 2 [a;(i7r + 6>)] sinh[7r(l - fe)x/2] 
xsinh7rxcosh(7rx/2) sinh(7r/ca;/2) 

1 



(55) 



d 



647r 3 Coo(i7r)C4(ivr) 
Thus the leading asymptotics of (I a I b ) and (I a I b ) are 

(r(r,x)/ fe (0,0)) = 2rf 2 5 afc x (56) 
J d^# 2 |F 2 (#i 2 )| 2 e (9l+e2) exp[-mr(cosh#i + cosh# 2 ) + imx(sinh#i + sinh# 2 )] 
(I a (T,x)l\0,0)) = 

S ab 2d 2 J d6 1 dd 2 \F 2 (6 12 ) | 2 exp [-mr (cosh Q x + cosh 6 2 ) + imx(sinh #1 + sinh 9 2 )] (57) 



where 



I^WI 2 - RW^iP^l IC " (fi|(il * (58) 



For the Fourier transform we have 

iD(s 2 \ = 

Vs 2 — 4m 2 



4r/ 2 

SmD(s 2 ) = - |^ 2 (g)| 2 , cosh 2 (#/2) = s 2 /4m 2 (59) 
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such that at the threshold we have the same behavior as for AT = 1: 

Q7B.D(s 2 ) ~ V s 2 — Am 2 



(60) 



Single-electron Green's function for the model of stripes 

All correlation functions discussed so far exist for both models depicted on Fig. 1. Now 
we are going to discuss the single-electron function which exists only for the model of stripes 
(Fig. lb). For AT = 1,2 we can calculate the spinon Green's function G and find the 
following expression for the spectral function: 

1 / m 2 \ 3/8 

A R (u, q) 2 7 , 2 J (61) 

u — vq \uj z — [vq) z — m z J 

and for the tunneling density of states: 

.cosh-V/m) d#cosh(3fl/8) 
pH ~ Jo (^/m-coshg)3/3 ( 62 ) 



Correlation functions of the staggered magnetizations 

As we have stated above, the correlation functions of staggered magnetizations factorize 



into a product of two kinds of correlation functions: G and D (see Eqs.fl33|,36)). The function 
D is essentially non-perturbative and therefore especially interesting. We have demonstrated 
that for AT = 1 this function is non-zero. For N = 2 with the periodic boundary conditions 
this function vanishes. This is related to the fact that for the periodic boundary conditions 
the singlet sector is critical. The operator expfi-v/^^] contains vertex operators <To; at 
criticality correlation functions of such operators with vertex operators of different chirality 
always vanish: (cr o"o) = 0. 

To understand this better it is convenient to write the staggered magnetizations in the 
appropriate basis. Thus we have 

N 1>s = a^Tr[a{g + g+)] ± i^ ] Tr[a(g - g + )} 

N 2i4 = 4 2) TA°(G + G + )] ± ifi 2) Ti[a(G - G + )} (63) 

where g, G are the SU(2) matrices. The fields do, fio remain critical. The fields g, G have to 
be decomposed into their vertex operators and the interaction involves sectors with different 
chirality. 
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From Eqs.(|63"D we have 

«NiN!» = (((Wo)) + ((IH>IM>))) ((Tr[^]Tr[^ + ]}) 

((NxNs)) = («a a » + (((Tr[^]Tr[^])) + <<Tr[^ + ]Tr[^ + ]>>) (64) 

where ((a a )) = ((^oVo)) ~ |(wr) 2 + x 2 | 1/8 . 

There are no reasons to think that correlation functions on the chains with the same 
parity vanish. Here solitons enter with Lorentz spin 3/16. The asymptotics of the correlation 
function is given by 



2 

p- 1/A J d9e 3e/8 e- mTCOshe+ixmsinhe = 2p~ 1/A J d9K 3/4 (2mpcosh9) (65) 
where p 2 = (vr) 2 + x 2 . The Fourier transform is (s 2 = u 2 + (vq) 2 ) 

4tt J dpp 3/4 J (sp) J dBK 3/4 (2mpcoshe) ~ 

J dfl(cosh fl)~ 7/4 F(5/4, 1/2, 1; -s 2 /4m 2 cosh 2 6) (66) 

After the analytic continuation iuo = uo + iO we get the imaginary part (now s 2 = uj 2 — (vq) 2 > 
Am 2 ) 

*>,,) ~ d.^v^.^ ^A V2, 1/4; 1 - , 2 ) (67) 

Close to the threshold s = 2m we have 

XV, 0) ~ (« 2 - 4m 2 )" 1 / 4 (68) 

At s 2 3> m 2 Eq.(^) gives the correct asymptotics x" ~ s_1 which indicates that the two- 
soliton approximation may give a reasonable description throughout the entire range of 
energies. Recall that for N — 1 the power was —1/2. It is tempting to speculate that the 
threshold singularity further diminishes with an increase of N. 

The case of four chains introduces some new features. Some of them, as we believe, are 
accidental and some are generic. The generic feature which persists for higher number of 
chains (see above) is the presence of singlet degrees of freedom. Here they appear in the form 
of the critical Ising model. The accidental features are the criticality of the latter model and 
the non-Abelian statistics of the massive kinks. These two properties are closely related and 
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are unstable with respect to a change in boundary conditions. In the most general setting, 
as we shall see in the next Section, the gap in the Ising model sector will be finite, but much 
smaller than the magnetic gap. 



V. TOWARDS INFINITE NUMBER OF CHAINS: A TRIADIC CLUSTER EX- 
PANSION. 



As the reader probably understands, the problem in question is difficult, naive attempts 
to develop a self-consistent expansion schemes fail (see Section IIIA) and we have to resort 
to some other methods. The exact results for two and four chains discussed in the previous 
Section give a glimpse of the complexity of the problem. Unfortunately, for higher number of 
chains there are no exact solutions. However, a non-perturbative analysis can be imagined, 
though for a somewhat modified model. In this modification the interactions are set up in 
such a way that the chains are assembled into clusters of three (triads), nine (enneads), 27 
chains etc. (see Fig.|| ). 
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FIG. 2: The hierarchical model for 36 chains. The circles with different color denote currents 
with different chirality. The green bonds are the strongest, the magenta ones are of intermediate 
strength and the yellow ones are the weakest. 
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The interaction inside of a cluster of a given size is supposed to be stronger than interac- 
tions between the clusters of larger sizes. We believe that the study of such modified model 
sheds light on classification of the excitations and establishes the structure of the effective 
action for the collective modes below the soliton mass gap. 



A. The fundamental triad (three chains) 

The solution for two chains with one chirality (let it be the left one) coupled by the 
current-current interaction to a chain with opposite chirality was described in |57| and |58 



(see also 59 



The Hamiltonian density is 

n 



W - A - , ™ TT 

-^XodxXo + y : II : 



),iU :JJ:+7lJ (69) 



3 

where I a and J a are right SU2(2) and left SUi(2) currents, respectively. The expression in 
the square brackets describes a sum of two chiral SUi(2) WZNW models. Such represen- 
tation was already discussed in Section IV.B. The right-moving Majorana fermion does not 
participate in the interaction. In the ultraviolet limit, the model represents a sum of three 
chiral conformal field theories: the k = 2 SU(2) WZNW model with the right central charge 
C R = 3/2, the right-moving free Majorana field with C R = 1/2 and the SUi(2) WZNW 
model with the left central charge Cl = 1- The current- current interaction generates a 
massless RG flow to the infrared critical point. The theory in the infrared is represented by 
two free Majorana fermions with opposite chirality xo an d Xo and the right-moving sector 
of the SUi(2) WZNW model: 

Hm = — : jj : +— (xodxXo ~ XodxXo) (70) 

Under the RG flow the fields transmute in the following way: 

I^2j + ..., ixoX ~> iXoXoj + ■•• (71) 

where the dots stand for less relevant operators. 



B. The first ennead (nine). 

Let us consider the case of nine chains with open boundary conditions, arranging them 
into three triads in such a way that the bare coupling constant inside each triad 70 is greater 
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+ 

C \ 
Ising model 

FIG. 3: The renormalization scheme for three chains. In the infrared limit the system is equivalent 
to a single chiral C=l chain and a critical (non-chiral) Ising model. The circles with different 
thickness denote currents with different chirality. 

than the bare coupling 71 between the triads (see Fig. 3). Then we can have a situation 
when the IR fixed point for a given triad is already achieved (the corresponding energy scale 
is A ~ J||7 exp[— 77/7]), well before the renormalized coupling for the triad-triad interaction 
becomes of the order of 1. 

According to Eq. flTT|) , the renormalized coupling between any neighbouring triads (let us 
denote them 1 and 2; recall also that they always have different chirality!) in the infrared is 

7i(A)jij 2 + 7i(A)(ixoXo)i(iXoXo)2jij2 + ••• (72) 

The first term reproduces original interaction fl7|). Its presence in the Hamiltonian also 
converts the second term into a marginal operator: 

7i(A)(ixoXo)i(iXoXo)2jij2 = 71(1X0X0)1(1X0X0)2 + 7i(A)(iXoXo)i(iXoXo)2 : J1J2 : 
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7i = 7i(A)(jJ 2 ) (73) 

Thus by integrating out the high energy degrees of freedom in the ennead of chains one 
generates at energies smaller than A (the new UV cut-off) the effective action which contains 
the original action for a triad of chains with a renormalized coupling constant 71 plus the 
action for three critical Ising models coupled by the products of energy density operators: 

3 ri I 

Wising = «(x3»X ~ Xd x x)r + 7i(xx)2[(xx)i + (xx)s] (74) 

r=l lZ J 

The Ising subsystem decouples from the magnetic one. The resulting model is similar to the 
0(3) GN model discussed above (see the discussion around Eq. (|47|) ) . The only difference is 
that the symmetry is broken down to U(l). The interaction can be written in terms of the 
currents 

7i[(xx)i(xx) 2 + (XX)2(XX) 3 ] ~ li_[I x I* + I y I y ] + l\\l z V (75) 

which makes the model somewhat similar to the anisotropic Thirring model ( I a are SU2(2) 
currents made of Xr's and the coupling constants r y± ~ 71 and 7|| = on the bare level). 
We believe (the details will be given in a separate publication), that this particular type of 
anisotropy vanishes in the strong coupling regime in the same way as it does in the C-sector 
of anisotropic Thirring model [^4|,||55||. In other words in the strong coupling regime the 
model ([74]) is equivalent to the 0(3) GN model. 

Thus on each step of the triadic real space RG we generate the 0(3) GN models. Their 
number is equal to the number of clusters on the given level. For example, if the initial 
number of chains is 3^ and the first GN model appears for a cluster of 9 chains, the number 
of copies of GN models on the first level is 3 N ~ 2 . Then the next levels contain S^ -3 , 3 N ~ 4 
etc. copies such that the total number of such models in both parity sectors is S^" 1 . Taking 
into account that each of them has 2-fold degenerate ground state, we get the ground state 
degeneracy 2 A ^ //3 , where Af = 3 N is the total number of chains. The resulting ground state 
entropy is three times smaller than the one which was obtained for the uniform case. 

Now we can discuss the thermodynamics. In our cluster expansion we have the following 
energy scales. First, there is a sequence of crossover scales A > Ai, ... > Ajv_i (magnetic 
gaps) corresponding to crossovers inside of each cluster. Then there are energy gaps of the 
interacting Ising models M 2 ,M 3 ,... which are formed by clusters of 9, 27, 81, etc. chains. 
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FIG. 4: A schematic picture of renormalization of nine chains. 



The first scale in this sequence is of the order of 



M2 ~ Ai exp[— ir/'-fi] 



(76) 



Since 71 ~ 7^ we deem these energy gaps to be much smaller than the magnetic gaps. We 
conjecture that this difference survives even in the limit when all interactions become equal. 

All these crossovers affect the temperature behavior of the specific heat. The first 
crossover occurs at temperature T ~ Ao which is the crossover temperature of a funda- 
mental triad. Above this temperature we have a bunch of non-interacting chiral Heisenberg 
chains each having the heat capacity linear in T: 

7rT 

C(T>A ) = SC, C = — (77) 

Of 

where S = 3 L\\ is the total area occupied by the system. At A > T > Ai we effectively 
have S^ -1 critical Ising chains with central charge C = 1/2 and the same amount of non- 
interacting chiral Heisenberg chains. As a result the slope of the specific heat drops by 
the factor of 2. When the entire sequence of magnetic crossovers is passed, that is at 
Aat_i > T > M2 the specific heat is given only by the Ising models: 

TFT 



„ iT,l 1 1 



I2v 



(78) 



that is remains the same as after the first crossover. Thus we can say that in the thermo- 
dynamic limit the singlets ocuppy a half of the original Hilbert space. 

When the temperature falls below the first singlet gap M2, one may say to the first 
approximation that the Ising modes with the largests gap cease to contribute. After a 
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FIG. 5: The hierarhy of interactions and the renormalization process for 27 chains. The strongest 
bonds are brown, the intermediate are gree and the weakest ones are yellow. The Ising variables 
are shown as pink rectangles. The red and blue circles of various sizes depict spin currents of 
different chirality. Their size increases with every step of the RG process. 

certain crossover the linear slope in the specific heat drops by the factor of 9. The specific 
heat remains linear in T until M 3 is reached where the slope falls by another factor of 3 etc. 
In our scheme where the distribution of coupling constants between the clusters of different 
size is rather arbitrary, it makes little sense to do more detailed calculations. The only thing 
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FIG. 6: A schematic picture of the temperature dependence of the specific heat 

we can say is that C ~ T above certain temperature and than experiences a fast decrease. 
In our cluster expansion there are two areas with linear specific heat characterised by the 
slopes differing by the factor of 2. This due to the fact that within this approach the gap 
for magnetic excitations A is different from the gaps for singlet excitations. Whether this 
feature will survive in the limit of uniformly coupled chains is open for debate. 

VI. SOME LEFTOVERS 

It would make a lot of sense to study the models described in this paper numerically. 
Since numerical calculations will probably be performed for systems with restricted number 
of chains, we decided to describe certain tractable cases in more detail. 

A. An approximate solution for N = 3 (six chains). 

In this case the magnetic subsystem is equivalent to AT = 1 (that is to the SU(2) Thirring 
model) and the Ising subsystem contains two Ising models with the interaction. The 
latter is equivalent to the spinless Thirring model with a truly marginal interaction. Thus 
the Ising sector is critical and is described by the C =1 Gaussian model. 
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B. An approximate solution for N = 6 (twelve chains). 



Applying the procedure outlined above to the system of twelve chains with periodic 
boundary conditions (see Fig. 3), we obtain the N = 2 model in the magnetic sector plus 



FIG. 7: A schematic picture of twelve chains. The fat lines correspond to larger exchange 
interactions. The circles of different size denote chains with different chirality. 

Thus the excitation spectrum includes GN 0(3) soliton mode, one gapless Majorana 
fermion and massive excitations from the Ising sector. 

VII. SINGLE-ELECTRON GREEN'S FUNCTION FOR THE MODEL OF 
STRIPES AND THE SPINON PROPAGATOR FOR N = oo. 

Now we shall try to use the wisdom accumulated in the exact solution of iV = 2 case 
to obtain results for the case of infinite number of chains. Though our ultimate goal is to 
calculate the correlation function of staggered magnetizations, the path to it lies through 
the single-electron Green's function, which exists only for the model of stripes (Fig. lb). 
The reason for this will become clear in the process of calculation. 

Let conduction electrons belong to even chains. Then the electron creation operator is 



where v c is the Fermi velocity. Though, in general, v c may be quite different from the spin 
velocity, we shall not consider this possibility. With v c = v s the model is (l+l)-dimensional 
Lorentz invariant. The charge fields ip c , (p c are free Gaussian fields and hence there are no 
correlations between the fields belonging to different chains. 






(79) 



Grr(t,x) 



1 



G 2n ,2n(r,x), G LL (r,x) 



1 



G 2 n,2n(r, -x) (80) 



(tv c + ix) 1 / 2 



(tv c — ix) 1 / 2 
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where function G was defined in Eg. (|32|) and can be called the spinon Green's function. 
Since the electron Green's function must be a single-valued function of x at r — ► and 
a single-valued function of r at x = 0, we conclude that the spinon function G must be 
a double- valued function to compensate for the double- valuedness of (rv c ± ix)~ 1 ^ 2 . This 
suggests that the spinons are semions, as it was suggested by Laughlin pOfl . 

Now we shall use these facts to determine the asymptotics of the single-fermion Green's 
function in a similar way as it was done for the (l+l)-dimensional Thirring model [fl}|,||44 



Namely, we shall combine these arguments with the Lehmann expansion for the Green's 
function. In this expansion we shall take into account only terms with an emission of a 
single massive soliton. The spectrum of this soliton is 

E(k, k ± ) = y/(vk) 2 + M 2 (k ± ) (81) 

where M(k±) ~ exp[— irv/'y} is as yet unknown function (the soliton gap). We believe that 
this gap always remains finite for any momentum. A convenient parametrization of the 
energy and the momentum component along the chain direction is 

E = M(k ± ) cosh 9, vk x = M(k±) sinh 9 (82) 

In these notations a Lorentz rotation on an 'angle' 7 corresponds to the shift of rapidity 9 — > 
# + 7. Since an operator with Lorentz spin S transforms under such Lorentz transformation 

as 

A s - e^ s A s 

We generalize for N = 00 the formulae obtained for N — 1 and N = 2 case treating the 
vanishing of D-function is the iV = 2 as an artefact of periodic boundary conditions. 

G ~ (vt - ix)~ 1/2 Z{k ± ) exp[-M{k ± )^T 2 + {x/v) 2 ] 



D ~ f^M^^K^M^r 2 + (x/v) 2 } (83) 

where Z, f are yet unknown functions. The /-function may vanish at some point in momen- 
tum space, as suggested by the example of four chains. Notice that both /c^-dependence of 
G and the very existence of D are effects exponential in I/7. These equations dictate the 
following asymptotic form for the correlation function of the staggered magnetizations: 

^±Z(k ± )Z(q ± + k ± ) exp{-[M(k ± ) + M(q ± + k ± )}p} (84) 
jAf(k ± )f(k ± + g ± )[M(^)M(g ± ] 1 /2^ [M(fc ± )p]^o[M(A; ± + q ± )p] 
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where p 2 = r 2 + (x/v) 2 and \i- is a correlation function between the chains with the same 
and X2 with different parities. 

VIII. CONCLUSIONS 

Let us make a summary of our results. 

• We have found models of magnets with a short range Heisenberg exchange which have 
neither magnetic nor spin-Peierls order. 

• We have given a formal prove that these models possess excitations with spin 1/2. 

• We established the existence of T = critical point and a massive ground state 
degeneracy. The ground state entropy is proportional to the number of chains. 

• We have also established that the low-energy Hamiltonian separates into two weakly 
interacting parts describing sectors with different parity. 

• We established that the correlation functions of staggered magnetizations can be writ- 
ten as real space products of correlation functions of non-local operators belonging to 
the sectors with different parity (see Eqs. (|32"|, |33|)). 

• We have solved the problem exactly for the case of two and four coupled chains (the 
latter one with periodic boundary conditions). The results are consistent with the 
statements made for the case of infinite number of chains. The spectrum contains a 
massive amount of singlet excitations. In the language of gauge theory these excita- 
tions describe dynamics of Wilson loops. 

A natural question is whether the models considered in our work can be realized. We 
believe that the the answer is positive. The model of stripes may well be relevant in a 
strongly underdoped regime of copper oxides. 

We have already mentioned in Introduction, that there are other models of fractionaliza- 
tion, such as suggested in 0,|§ together with the corresponding dimer models mentioned 
earlier. Though these works use Hamiltonians not based on any microscopic electron models, 
one might hope that they capture some features of the solution. In particular, it would be 
interesting to study in detail how singlet excitations, which in our model are associated with 
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the 0(3) GN models emerging on each triad of chains, are related to Z2 vortices (visons) 
introduced in pH , [ 6"I] (see also the recent work ||62|| , which essentially clarifies the concept of 
vison). We leave this question, as many others, to the future research. 
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X. APPENDIX A. GAUGE THEORY FORMULATION. 

The relationship with the gauge theories can be inferred from the fact that the level 
k — 1 SU(2) Kac-Moody currents can be identified with (iso)spin currents of massless Dirac 
fermions: 

a a being Pauli matrices. The Hamiltonian (^),(0) becomes 

7 

As gauge theories of the RVB state, this theory possesses redundant charge degrees of 
freedom which do not participate in the interactions. Since every approximation violates 
this subtle property, dealing with the charge sector would become the same awkward problem 
as it is in the standard approach to the RVB gauge theory, once one decides to adopt this 
fermionic representation. Let us, however, follow the well-trodden path for a while just to 
make sure that the model we are discussing does fall in the category of RVB liquids. To this 
end, we use the identity 

o$o% = 2P 12 -l, 
34 



where Pu is the permutation operator, and apply the Hubbard-Stratonovich transformation 
to rewrite the interaction term in Eq.(|86|) as 

\^fi+n,n( x ) | 2 



+ ... (87) 



27 

where the dots stand for the terms we deem irrelevant. The procedure essentially coincides 
with the conventional decoupling scheme in the RVB approach. The fact that such decou- 
pling here is done only in one lattice direction is not important provided one can justify that 
fluctuations of |A| may be neglected. As we shall see, excitations associated with breaking 
of singlets carry the largest spectral gap. This justifies the assumption about small fluc- 
tuations of |A|. Once the amplitude |A| ~ Jn exp(— nv/'-f) is frozen, one is left with the 
compact U(l) lattice gauge theory in the strong coupling limit (indeed, the gauge field has 
no bare F^-term which corresponds to infinite bare charge.). The vector potential A y is 
represented by the phase of A n ^ n+fJi . Since by omitting the term with R+R a L+L a we violate 
the decoupling of the charge degrees of freedom, we have to enforce the constraint on the 
absence of charge fluctuations by introducing the time component A of the gauge field. In 
the mean field approximation (that is, when fluctuations of the gauge fields are neglected) 
we obtain the spectrum of the 7r-flux state [ 12" ]: 



E(p) = yj(vk\\) 2 + 4|A| 2 cos 2 k ± (88) 

which has two Dirac-like conical singularities in the Brillouin zone. From what we have done 
in this paper it appears very dubious that this result will survive inclusion of fluctuations. 



XI. APPENDIX B. THE STRUCTURE OF STRONG COUPLING REGIME AND 
THE GROUND STATE DEGENERACY. 

Let us assume that the system is two-dimensional and the number of chains is even, 2N. 
We will drop the Lorentz-noninvariant part, C 1 ™ 1 , in which case one has a decomposition 

>~"n '-'n T -<- n • 

Let us first consider J2n^n- m the strong-coupling ground state, the following combina- 
tions of the chiral fields are locked: 

fi = ±l, m { £ =0,±1,±2,... (89) 



V 9 2n + ^P2n+n 



1 + 2mg 
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Assuming that periodic boundary conditions are imposed in the transverse direction, the 
set of equations (B9|) can be rewritten as 



V?0 + &2N-1 

¥2 + fil 

^2 + <f3 

<P2N-2 + l ^2N-Z 
<P2N-2 + <^27V-1 



7T r 



1 + 2ml 



(-) 



1 + 2m, 



(+) 


(-) 



1 + 2m. 



1 + 2m 2 +) 



71 r 
8 

W r 



1 + 2rri2]J-2 
1 + 2m2jv_2 



(90) 



Considering pairs of neighboring equations in (^) one finds that the integers m 2 ^ satisfy 
the condition: 

N-l N-l 

(91) 



N-l N-l 
\- (+) \ - (-) 
m 2k = 2^ ™>2k ■ 



k=0 k=0 

With this constraint, the number of independent fields, locked in the ground state, is 2N—1. 
To select these fields we first define N scalar fields 



Xn = </>2n + <^2n+l, 71 = 0, 1, 2, N - 1 



(92) 



together with their iV dual counterparts 



$n = (f2n ~ ¥2n 



According to (0), the fields \n are locked at 



+1- 



l + 2m« 



(93) 



(94) 



On the other hand, there are N more combinations, (f2 n + y?2n-i ; that get frozen, again 
according to ( [39]) . We can express them in terms of \ n and d n : 



<f2n + f2n-l = ~ (Xn + $n + Xn-1 ~ #n-l) 



l + 2m 



(-) 

2/i 



(95) 



Thus iV relative dual fields 



H = _(0 n -0 n+ i), n = 0,l,2,...,iV-l 



(96) 
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are also locked: 



Q(-) 



7T r 



.(+) i ™(+) 



,(-) 



m 2n + m 2n'+2 ^ m 2n+2 



(97) 



Notice that transverse periodic boundary conditions imply that 



(which is actually the same condition as Eg. (pi]) ) and so the number of independent relative 
fields is N — 1. We choose them to be 



U , Vi , ... V 



N-2 



(98) 



The total field 



JV-l 



9* 



71=0 



(99) 



remains unlocked and, hence, disordered. 



Thus the N fields Xm Eq.( ^) , are locked while N their dual counterparts $ n are not; 
only thier N —1 combinations (|98|) are. The frozen values of these 2N — 1 independent fields 
characterize the vacuum state of J2 n ^n- 

The same can be done for J2 n £-n : we introduce 



■0n = V?2n + f2n+li 



(100) 



and select iV fields ip n and iV — 1 relative duals 



..(-) ..(-) ..(-) 

U! , U! 1 , ... , UJ N _ 2 , 



(101) 



where 



The total field 



(-) _ ~ ^n+l 

V2 ' 



JV-1 

E 

n=0 



remains disordered. The ip n — fields are frozen at 



7r r 



1 + 2m ( J +1 



(102) 
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while ui > - at 



uJ 



{-) 



7T r 



2m 



(+) 

2n+l 



m 



(-) 

2n+l 



m 



H 1 

2n+3 



(103) 



We can now express the physical fields of individual chains in terms of the locked fields 
and two extra fields that remain disordered. We have: 



$2n = 


n (Xn + + - W n ) , 


&2n = 


2 (Xn - + $n + LU n ) , 


$2n+l = 


^ (Xn + Ipn ~ fin + U n ) , 


@2n+l = 


2" (~Xn + + + 



(104) 

The elements of the Wess-Zumino matrix field g n (x) contain exponents e ±n/ ^* n and 
e ±l ^" en . According to Eqs. ( |104j ), in the strong coupling phase these exponents are pro- 
portional to 

ew { i \I¥ e *) and exp (Vl^*)' 

respectively, and thus will have vanishing expectation values. This effect, however, disap- 
pears in the thermodynamic limit. 



Transverse dimerization and degeneracy of the ground state 

Spontaneous transverse dimerization has been identified for the case of two chains in 
Ref . ||25|| . In the continuum limit, the corresponding order parameter is given by a simple 
expression: N n (x) • N n+1 (x). For the 2N-chain model 



«2» • -^2» + J ~ COS V27r(6 2 n — @2n+l 
1 



+ 2 



COS V27r($ 2n - $2n+l) - COS V27r($2n + $2n+l) 

1 r 



cos V27r(xn - VvJ + 2 [ cos v27r(i? n - u n ) - cos V27r(x n + ip n ) 



Using Eqs. (|94]) and (|102[) when averaging over the ground state, one finds that 



??2n = (N 2n -N 2n+ l) ~ (cOSV27r(Xn-^n)) --(cOSV / 27r(Xn + V ; n)) 

OC ^(-1)^- "W 



(105) 



38 



We observe the existence of two, doubly degenerate values of this local order parameter; its 
sign depends on the parity of the integer m^ n — . 
Similarly 



N 2 „+l ■ N 2n+2 ~ COS V / 27T(02n+l - @2n+2 
1 



+ 



/7T 

~ COS t/~ 



COS V27T($2n+l - $2n+2) ~ COS V 27r($ 2 n+l 
n ~r Xn+1 



2n+2j 



1 [W 



1 

- COS . , 

2 V 2 



On + VWl) + (Xn + Xn+l) ~ (0 n - CJ n ) 



(106) 



The first and third terms in the r.h.s. of (|106|) have nonzero expectation values. Using 
the locked values of the corresponding fields given above, we find that arguments of these 



two cosines are ir 



V2n+1 = (N 2n +1 • N 



(+) 
2n+l 



+ m 



(-) 

2n+2 



+ 1 



respectively. Therefore 



L>; i+ 2) OC ^(-l) m 2"+l 



(107) 



Inspecting ( |105| ) and ( |107| ) we observe that the two signs of local order parameters r] m 
reflect a spontaneously broken Z 2 symmetry. This is the symmetry related to independent 
translations by one lattice spacing along the chains. Notive, however, that for different pairs 
of chains the signs of the order parameters are uncorr elated. The only condition imposed 
on the 77's is 

2N 

-l) Q = 1 (108) 



2N 
n=0 



where the integer 



iV-l 



2m+l\ 



vanishes according to relations ( |9*ID in the (+) sector and their counterparts in the (-) 
sector. Eq.( |108j ) is not a restrictive condition; it simply says that the numbers of positive 
and negative 77's should be even. Thus the ground state of the system exhibits a huge degree 



of degeneracy. Taking the constraint (|108|) into account, the number of the degenerate 
ground-state configurations in a system of 2N chains is estimated as 



2N 
n=0 



2iV 



2iV 

E 

n=0 



(2N)\ 



i2N-l 



n\{2N -n)\ 
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